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ABSTRACT 

Aims. We performed deep optical observations of the field of an old, fast-moving radio pulsar PSR Bl 133+16 in an attempt to detect 
its optical counterpart and a bow shock nebula. 

Methods. The observations were carried out using the direct imaging mode of FORS1 at the ESO VLT/UT1 telescope in the B, R, and 
H a bands. We also used archival images of the same field obtained with the VLT in the B band and with the Chandra/ACIS in X-rays. 
Results. In the B band we detected a faint (fl=28?l±0?3) source that may be the optical counterpart of PSR Bl 133+16, as it is 
positionally consistent with the radio pulsar and with the X-ray counterpart candidate published earlier. Its upper limit in the R band 
implies a color index B-R^U?5, which is compatible with the index values for most pulsars identified in the optical range. The derived 
optical luminosity and its ratio to the X-ray luminosity of the candidate are consistent with expected values derived from a sample of 
pulsars detected in both spectral domains. No Balmer bow shock was detected, implying a low density of ambient matter around the 
pulsar. However, in the X-ray and H„ images we found the signature of a trail extending ~ 4" - 5" behind the pulsar and coinciding 
with the direction of its proper motion. If confirmed by deeper studies, this is the first time such a trail has been seen in the optical 
and X-ray wavelengths. 

Conclusions. Further observations at later epochs are necessary to confirm the identification of the pulsar by the candidate's proper 
motion measurements. 

Key words, pulsars: general - pulsars, individual: PSR Bl 133+16 - stars: neutron 



^ ■ 1. Introduction The non-thermal component, which also is frequently ob- 

served in a wider spectral range including X-rays, is believed to 



Until now optical emission has only been detected from < 1 % be powered by the NS rotational energy loss E, called a spin- 
el 'of ^1500 known radio pulsars (e.g., Mignani et al. 2005). down luminosity. The parameter 77 = L/E, where L is the radia- 
Even such a small number of optical identifications implies that tive luminosity, describes the efficiency of the transformation of 
rotation-powered neutron stars (NSs) can be active in the opti- the rotational energy into the emission. The most striking result 
cal, as well as in the radio range. This follows from the fact that of the optical study of ordinary pulsars is that the two old pulsars, 
a small group of the optical pulsars contains not only young, PSR B1929+10 and PSR B0950+08, with rather low E, have al- 
~ 1 kyr, and energetic objects like the Crab pulsar (Percival et most the same optical efficiency as young and energetic Crab- 
al ll993l >, but also much olde r pulsars such as -3 .1 Myr PSR like o bjects with much higher spin-down luminosities (Zharikov 
B1929+10 (Pavlov et al. 1 19961 Mig nani et al.[2002 ]>, -17 4 My r et al. 120041 At the same time, the efficiencies of middle-aged 
PSR B0950+08 (Pavlov et al. 1 19961 Zharikov et al. 120021 12004). pulsars are significantly lower. In addition, a strong correlation 
and even the very old ~1 Gyr recycle d millisecond pulsar between the non-thermal optical and X-ray luminosities of pul- 
PSR J0437-4715 (Kargaltsev et al. 120041) . The power-law non- sars was found, indicating a general o rigin of the emission in 
thermal spectral component is dominant in the pulsar's optical both ranges. (Zharikov et al. 2004 2006). To confirm these find- 
emission, presumably originating in the magnetospheres of the ings it is necessary to increase the number of optically identified 
NSs. pulsars. 

An old, ~5 Myr, nearby pulsar PSR Bl 133+16 has al- 



Send offprint requests to: S. Zharikov, most the same parameters as the two old objects mentioned 

e-mail: zhar@astrosen.unam.mx above (see Table Q). It is located at a high galactic latitude, 

* Based on observations made with ESO telescope at the Paranal I = 242°, b = 69°, implyi ng a l ow interstellar extinction 

Observatory under Program 7 1 .D-0499(A) and with archival ESO VLT E(B - V) = 0.04 (Schlegel et al. |1998) . The direct proper motion 

data, obtained from the ESO/ST-ECF Science Archive Facility. and annual parallax measurements in the radio range by Brisken 



2 S . Zharikov, et al. : Optical observations of the PSR B 1 1 3 3 + 1 6 



Table 1. Parameters of several nearby old pulsars. 



Pulsar 


age 


distance* 


parallax* 


AC 




v*± 


log(£) 


optical band 


magnitude 




log (t ) 


pc 


mas 


mas yr 1 


mas yr 


km s 


ergs s 






PSR B 1929+ 10 


6.49 


331(10) 


3.02(9) 


94.8 


43 


163(5) 


33.59 


F342W 


26.8° 


PSRB1133+16 


6.70 


350(20) 


2.80(16) 


-73.95 


368.05 


631(38) 


31.94 




2S.l d 


PSR B0950+08 


7.24 


262(5) 


3.82(7) 


-2.09 


29.5 


36.6(7) 


32.74 


B 


27.07* 


PSR JO 108- 1431 


8.23 


60-130 




< 26 


<78 


< 50bo 


30.75 


B 


> 28.6'' 



Brisken et al. 2002, " Pavlov et al. 1996; * Zharikov et al. 2002, 2004; c Mignani et al. 2003 ; d this paper 



Table 2. Log of the VLT/FORS1 observations of the PSR B1133+16 field during the 2003-2004 periods 



Date 


UT 


Band, exp. 


Exp. 


secz 


Seeing 


Sky 


Date 


UT 


Band, exp. 


Exp. 


secz 


Seeing 


Sky 


UT 




number 


sec. 




arcsec. 


counts 
sec pix 


UT 




number 


sec. 




arcsec. 


counts 
sec pix 


03 Apr 2003 


02:59 


B 1 


1020 


1.33 


0.95 


3.04 




08:03 


B 12 


1020 


1.57 


0.97 


2.53 




03:17 


B 2 


1020 


1.31 


0.92 


3.02 


25 Dec 2003 


07:56 


R7 


485 


1.58 


0.57 


17.38 




03:35 


B 3 


1020 


1.31 


0.91 


2.98 




08:05 


R8 


485 


1.54 


0.55 


16.80 


J> 


03:55 


R 1 


750 


1.32 


0.77 


27.24 




08:15 


B 13 


30 


1.50 


1.01 






04:08 


R2 


750 


1.33 


0.71 


25.78 


19 Jan 2004 


07:24 


H„ 1 


1550 


1.37 


0.67 


0.65 




04:22 


R3 


750 


1.34 


0.71 


24.75 




07:51 


H„2 


1550 


1.33 


0.64 


0.67 




04:35 


R4 


750 


1.37 


0.73 


24.21 


25 Jan 2004 


07:08 


H„ 3 


1550 


1.35 


0.83 


0.71 


04 Apr 2003 


03:33 


R5 


750 


1.31 


0.91 


17.16 




07:35 


H„4 


1550 


1.32 


0.89 


0.74 




03:46 


R6 


750 


1.32 


0.84 


17.47 


26 Jan 2004 


05.46 


H„5 


1550 


1.61 


0.93 


0.59 


09 Apr 2003 


02:37 


B 4 


1020 


1.32 


0.87 


5.24 




06:13 


H„6 


1550 


1.48 


0.83 


0.57 




02:55 


B 5 


1020 


1.31 


0.86 


4.14 




06:41 


H„l 


1550 


1.40 


0.99 


0.59 




03:12 


B 6 


1020 


1.31 


0.95 


3.08 




07:08 


H„ 8 


1550 


1.35 


1.11 


0.57 




03:32 


B 7 


1020 


1.32 


0.95 


2.65 


27 Jan 2004 


07:22 


H ff 9 


1550 


1.32 


0.88 


0.75 




03:49 


B 8 


1020 


1.33 


0.81 


2.62 


31 Jan 2004 


07:03 


R9 


750 


1.33 


0.56 


23.42 




04:08 


B9 


1020 


1.36 


0.79 


2.70 




07:17 


R 10 


750 


1.32 


0.66 


15.52 


24 Dec 2003 


07:28 


B 10 


1020 


1.79 


0.94 


2.56 




07.30 


R 11 


750 


1.31 


0.89 


15.15 




07:46 


B 11 


1020 


1.67 


0.88 


2.52 




07:43 


R 12 


750 


1.31 


0.83 


14.95 



et al. (120021 1 yield a high transverse velocity of 631+30 km/s and 
a short distance to the pulsar 350+20 pc. The pulsar is younger 
than PSR B0950+08, but its spin-down luminosity, 8.8 x 10 31 
ergs s -1 , is about an order of magnitude lower than those of PSR 
B 1929+20 and PSR B0950+08. Nevertheless, it is higher than 
that of the nearby old PSR J0108-1431, whose optical counter- 
part has been unsuccessfully searched by Mignani et al. (2003). 
The high transverse velocity is promising for detecting an H a 
bow shock nebula expected to be produced by the supersonic 
motion of the pulsar in the interstellar matter, as has been found 
around several rapidly moving pulsars and radio-silent NSs (e.g., 
Gaensler & Slane 2006). Recently, the field of PSR Bl 133+16 
has been observed in X-rays with the Chandra/ACIS and a 
faint X-ray counterpart candidate was found at a flux level of 
(8+2)xl0~ 15 erg s^'cirT 2 in 0.5 - 8 keV range (Kargaltsev et 
al. [2006). Given this value and using an empirical relation be- 
tween the optical and X-ray luminosities of pulsars (Zharikov et 
al. 2004 2006; see also Zavlin & Pavlov 120041) . one can expect 
to detect the pulsar in the optical range at a sensitivity level of 
27-29 magnitudes. 



In this paper we present the results of a deep optical imag- 
ing of the PSR Bl 133+16 field taken with the ESO Very Large 
Telescope (VLT) in B, R c , and H a bands to search for the opti- 
cal counterpart and the bow shock nebula of the pulsar. We also 
used archival VLT images of the field obtained earlier in the B 
band and the X-ray data from the Chandra archive. The obser- 
vations and data reduction are described in Sect. [2] Astrometric 
and photometric referencing are given in Sect. [3] We present our 
results in Sect. [4] and discuss them in Sect. [5] 



2. Observations and data reduction 

The observations were carried out with the FOcal Reducer and 
low-dispersion Spectrograph (FORSlQ) at the UT1 (ANTU) unit 
of the ESO/VLT during several service-mode observational runs 
from the beginning of April 2003 to the end of January 2004. 
A standard resolution mode was used with an image scale of 
~072/pixel and field of view (FOV) ~6:8x6.'8. The total inte- 
gration time was 12240 s in B, 8470 s in R c , and 13950 s in 
H a bands. The log of the observations is given in Table [2] To 
complement our analysis, we retrieved from the ESO ar chivefl 
VLT B-band images of the same field obtained under Program 
66.D-0069A of Gallant et al. (120001 in January 2001. These 
unpublished data were taken with the same telescope unit and 
instrument in a high resolution mode with an image scale of 
~ 071/pixel and ~3'4x3.'4 FOV. The total integration time was 
about 15 ks. The observational log is given in Table[3] Generally, 
all the observations were made under photometric conditions, 
instrumental magnitude variations of field stars through dataset 
being insignificant, but the sky background level, seeing, and at- 
mospheric extinction differ. For instance, the individual images 
in the B and R c bands for the 2003-2004 period can be subdi- 
vided into two distinct groups with a high and low background 
level (Table The H a data are more uniform in this respect, 
but seeing conditions vary from 0767 to l'.'ll. The 2001 data 
were obtained under better and less variable seeing conditions, 
but the sky background level was high for the exposures taken on 
January 18 and 19, 2001 during moonlight time. (Table[3) The 
data on January 25, 2001 were obtained under excellent seeing 
and sky background conditions, but the atmospheric extinction 

1 For the instrument details see www.eso.org/instruments/fors/ 

2 archive.eso.org 
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Table 3. The same as in Table [2]but for the 2001 period 



Date 


UT 


Band, exp. 


Exp. 


secz 


Seeing 


Sky 


Date 


UT 


Band, exp. 


Exp. 


secz 


Seeing 


Sky 


UT 




number 


sec. 




arcsec. 


sec pix 


UT 




number 


sec. 




arcsec. 


sec pix 


18 Jan 2001 


07:41 


B 1 


1000 


1.34 


0.51 


1.63 


-"- 


07:58 


B 10 


1000 


1.33 


0.83 


0.51 




07:59 


B 2 


1000 


1.32 


0.57 


1.63 




08:16 


B 11 


1000 


1.32 


0.92 


0.53 




08:16 


B3 


1000 


1.31 


0.64 


1.66 


25 Jan 2001 


06:56 


B 12 


1000 


1.37 


0.58 


0.61 




08:38 


B4 


2 


1.32 


0.65 


6.5 




07:14 


B 13 


1000 


1.34 


0.63 


0.60 




08:39 


B5 


20 


1.33 


0.65 


6.9 




07:31 


B 14 


1000 


1.32 


0.67 


0.59 


19 Jan 2001 


08:09 


B 6 


1000 


1.32 


0.79 


0.96 




07:50 


B 15 


1000 


1.31 


0.79 


0.58 




08:27 


B 7 


1000 


1.31 


0.85 


0.97 




08:10 


B 16 


1000 


1.32 


0.66 


0.57 




08:44 


B 8 


1000 


1.31 


0.85 


1.30 




08:28 


B 17 


1000 


1.33 


0.66 


0.57 


22 Jan 2001 


07:41 


B9 


1000 


1.32 


0.81 


0.51 

















Table 4. The list of USNO-B1.0 stars used for astrometrical ref- 
erencing with coordinates (epoch J2000.0) and their errors. 



Name RA DEC <x RA <r DEC 

hh mm ss dd mm ss mas mas 



was by ~0.16 magnitude higher than the standard VLT value, 
which may be caused by some cirri. 

The data reduction including bias subtraction, fiat-fielding, 
removing cosmic ray traces and alignment of each individual 
image to a reference frame was performed for all the data using 
standard IRAF and MIDAS tools. Taking the data non-uniformity 
into account, we considered different combinations of individual 
exposures in each band to get co-added images of the best qual- 
ity, deepness, and spatial resolution. A simple sum of all avail- 
able images for a given band and period was finally chosen as 
optimal and subsequently analyzed. Only very short, % 30 sec, 
exposures in the B band (see Tables[2]and[3]l were excluded. As a 
result, the FWHM of a stellar object in the composed images was 
«0775 and «1" in the B band for the 2001 and 2003-2004 pe- 
riods, respectively. The respective values in the R c and H a bands 
were »:0785 and The panoramic view of the pulsar field as 
seen with the VLT in the B band is shown in Fig.Q] 



3. Astrometry and photometry 

3.1. Astrometry 

The B-band image obtained by co-adding all B band frames 
taken during the 2003-2004 period was chosen as a base image 
for our astrometry because its FOV is twice bigger than in the 
2001 period. To compute a precise astrometric image solution, 
the positions of objects selected from the USNO-B 1 astrometric 
catalogP] were used as a reference. There are about 50 USNO-B 1 
objects in the FOV in contrast to none from UCAC2. The recent 



3 USNO-B 1 is currently incorporated into the Naval Observatory 
Merged Astrometric Dataset (NOMAD), which combines astrometric 
and photometric information of Hipparcos, Tycho-2, UCAC, Yellow- 
Blue6, USNO-B, and the 2MASS, www.nofs.navy.mil/data/fchpix/ 



release of the Guide Star Catalog (GSC-II v2.3.2)Q contains a 
similar number of stars in this FOV but provides no information 
on proper motions, and the declared astrometric errors (073) are 
higher than nominal 072 uncertainty of USNO-B 1 . We discarded 
the reference stars with significant proper motions and catalog- 
positional uncertainties k 074 along with saturated ones. The 
pixel coordinates of 30 objects considered to be suitable astro- 
metric reference points were computed making use of the IRAF 
task imcenter. The IRAF tasks ccmap/cctran were applied for 
the astrometric transformation of the image. Formal rms uncer- 
tainties of the initial astrometric fit were rather large, A RA = 
07287 and A DEC=07251, as compared with the 072 pixel scale 
of the image. Using standards with the smallest catalog uncer- 
tainties and fit residuals significantly improved the fit. For nine 
reference stars listed in Table [4] and marked in Fig. [TJ the for- 
mal rms are A RA = 07051 and A DEC=07097, which is smaller 
than the pixel scale. As can be derived from the Table|4] the local 
USNO-B 1 catalog rms uncertainties of the list in RA and DEC 
are «07045 and «07039, respectively, which is a factor of four 
smaller than the nominal catalog uncertainty. To refine the fit, 
we further removed standards with the largest catalog uncertain- 
ties, i.e. the first and the fourth ones from Table|4] This provided 
much smaller rms, A RA = 070278 and A DEC=07068, and in- 
dividual standard residuals ;S 0709, which is consistent within 
2-cr with the local rms of the catalog and with the positional un- 
certainties of the standards in the image (^ 0705). We adopted 
this fit as a final one. It is important that by starting with 13 stan- 
dards the image transformation became very stable and practi- 
cally independent of removing other standards with the largest 
residuals. At these steps coordinate reference points clustered 
together only within the same image pixel, gradually migrating 
to some limiting position within this pixel. The same result was 
obtained using the list of stars with coordinates from GSC-II 
v2.3.2. Therefore, we are quite confident the derived coordinate 
reference that is reliable at least at the level of the nominal mean 
catalog accuracy of «072. Combining that with the best fit rms, 
we obtained secure 1-cr uncertainties of our astrometric refer- 
encing in the RA and DEC as 07202 and 07211, respectively. 
They are comparable with the pixel scale and much smaller than 
the seeing value. 

The co-added 2001 B-band image was rebinned by 2x2 pix- 
els to get the same pixel scale as for the 2003-2004 period. All 
co-added 2001 B, R c , and H a images were aligned to the base 
one using a set of suitable unsaturated stars with an accuracy of 
better than 0.025 pixel size, and the above coordinate reference 
of the base image was adopted for these images. Formal errors 
of this referencing are negligible (^0'.'005) in comparison with 
the astrometric referencing uncertainties of the base image. The 



www-gsss.stsci.edu/Catalogs/GSC/GSC2/GSC2.htm 



1058-0205218 
1058-0205222 
1058-0205225 
1058-0205226 
1058-0205233 
1058-0205234 
1058-0205238 

1058- 0205246 

1059- 0203095 



11 35 52.8980 
11 35 54.9493 
11 35 58.2340 
11 35 58.7207 
11 36 01.6207 
11 36 01.9007 
11 36 03.9327 
11 36 10.1713 
11 36 14.1813 



+ 15 52 07.000 
+ 15 52 03.040 
+ 15 51 19.070 
+ 15 49 47.110 
+ 15 49 57.610 
+ 15 49 40.620 
+ 15 52 32.970 
+ 15 53 21.120 
+ 15 54 37.130 



192 


144 


25 


23 


125 


239 


307 


138 


98 


97 


46 


47 


34 


54 


4 


19 


74 


109 



4 
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Fig.l. PSR Bl 133+16 field as seen in the B band with the VLT/FORS1 using a standard resolution mode during the 2003-2004 
period. The position of the pulsar is marked by +. The large box marks the FOV of high-resolution mode observations during the 
2001 period. The x-shaped crosses and numbers mark the USNO-B1.0 reference stars from Tableland secondary photometric 
standard stars from Table[7] respectively. The region within a small box is enlarged in Figs. [2] and [4] 



selection of another summed image, R c or H a , as a base does not 
change the result. 

The radio positions of PSR Bl 133+16 at the epochs of the 
VLT and Chandra observations (Table|5]l were determined using 
radio measurements of the pulsar proper motion made with the 
VLA by Brisken et al. 120021 The position errors in Table [5] ctra 
and ctdec, include the errors of the radio position at the reference 
epoch J2000 (15 mas for both coordinates) and the uncertainties 
in the pulsar proper motion (0.38 mas and 0.28 mas in RA and 
Dec, respectively). 

Combining the errors of our astrometry and the radio posi- 
tion errors, we derived uncertainties of the pulsar position in our 
images. Observations in different bands were splinted by differ- 
ent sets on a time basis from a week to 8-9 months, and the pul- 
sar has a large proper motion, particularly along the declination 



(Table Q]) (cf. Tables [2] and [3), therefore one has to account for 
a systematic uncertainty related to a shift in the pulsar position 
during this time. It is mostly significant for the B and R band 
observations of the 2003-2004 period when the pulsar shifts by 
~0'.'197 and ~0'.'229, or by about 1 CCD pixel, from the begin- 
ning to the end of the observations. We added these shifts as 
uncertainties squared and calculated the expected optical coun- 
terpart coordinates as a simple mean of those at the beginning 
and the end of the observations for given resulting frames and 
periods. The results are summarized in Table [6] As seen, de- 
spite a rather accurate astrometric referencing of the summed 
images, the resulting 3cr position uncertainty of the expected op- 
tical counterpart can be as large as 1 arcsecond in the DEC due 
to the large proper motion of the pulsar. 
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Table 5. PSR Bl 133+16 radio coordinates at the epochs of the VLT and Chandra observations. 



5 



Epoch 


RA 


DEC 


O-RA 


<Tdec 


2450000+ 


hh mm ss 


dd mm ss 


mas 


mas 


1544.5/2000* 


11 36 03.1829 


+ 15 51 09.7257 


15 


15 


1932.5/2001 Jan. 


11 36 03.1774 


+ 15 51 10.1189 


15 


15 


2738.5/2003 Apr 


11 36 03.1662 


+ 15 51 10.9228 


15 


15 


2997.5/2003 Dec. 


11 36 03.1625 


+ 15 51 11.1898 


15 


15 


3029.5/2004 Jan. 


11 36 03.1621 


+ 15 51 11.2221 


16 


15 


3424.5/2005 Feb. 


11 36 03.1565 


+ 15 51 11.6201 


16 


16 


* Brisken et al. 2002 



Table 6. Expected optical counterpart coordinates at the 
summed VLT frames. 



Frame/ 


RA 


DEC 




3cr dec 


Period 


hh mm ss 


dd mm ss 


arcsec 


arcsec 


S/2001 


11 36 03.1774 


+ 15 51 10.1189 


0.61 


0.63 


B/2003 


11 36 03.1644 


+ 15 51 11.0563 


0.63 


1.02 


fl,../2003-04 


11 36 03.1642 


+ 15 51 11.0719 


0.63 


1.10 


/4/2004 


11 36 03.1621 


+ 15 51 11.221 


0.61 


0.63 



3.2. Photometric calibration of the B and R c images 

For precise photometric calibration of all datasets taken in differ- 
ent nights and periods, a two-step approach known as "differen- 
tial photometry" was applied. At first, absolute magnitudes of a 
set of relatively bright stars visible in all target frames (so called 
"secondary standards") were derived accurately. This was done 
using images of the primary Landolt's standards obtained during 
the same night chosen as a "reference". In the second step the co- 
added target images were directly calibrated using the secondary 
standards. An advantage of this approach is that it is not neces- 
sary to consider possible magnitude zero-points and/or extinc- 
tion variations from night to night since the secondary standards 
and the target are in the same frames. 

The night of 02/03 April 2003 was used as a reference for 
photometric calibration and, for this night, the zero-points^ in 
the Johnson-Cousins B and R c bands, and color terms establish- 
ing the relation between the real B and R c magnitudes and re- 
spective instrumental magnitudes mg and m R were derived using 
1-43 Landolt's standar d stars from Rubinl49, PG1047+003, and 
SA1 10 fields (Landolt [T992j l: 

B = 27.17(3) + m B -0.04(2)(m B -mR) (1) 
R = 27.38(2) + m R + 0.01(2)(m B - m R ). 

For a star-like object, the instrumental magnitude in j-th band 
is defined as 

mj = -2.5 log{fi p lt exp ) - 6m j - kj secz , (2) 

where f„ p is the source flux in counts for a given aperture, t exp the 
exposure time, kj the extinction factor, secz is airmass, and dntj 
the correction for a finite aperture derived from a point spread 
function (PSF) of bright stars in the image. The signal-to-noise 



5 The zero-points were calculated for the flux in ADU/s units. The 
Paranal Observatory web page lists this night as a photometric one and 
our zero-points, when converted to electron/s units, are consistent with 
those provided for this night by the VLT/FORS 1 team. 

6 A maximum number of standards per night observed during our 
program. 



ratio S /N and magnitude uncertainty Am are calculated in a stan- 
dard way as (Newberry |1991| l 



N ^fap/g + n ap T 2 bg (l + l/n bg ) 

Am = 1.0856 (^j , (4) 

where crb g is the standard deviation of the background in counts, 
n ap the number of pixels in the source aperture, g the gain, and 
rib g the number of pixels in area used for the background mea- 
surement. With our seeing conditions, the source aperture with 
a radius of 3 pixels (at the image scale wO'.'2/pixel), for which 
S IN reaches a maximum, was found to be optimal for stellar ob- 
jects in all our images. The atmospheric transparency during our 
reference night was close to a standard value and we used the 
nominal extinction factors kg = 0.25 and k R = 0.08 provided for 
the VLT sitefl 



Table 7. Secondary standard stars used for photometric refer- 
encing. 



NN 


B 


R c 


Hi" 


NN 


B 


R c 


Uf 




mag 


mag 


mag 




mag 


mag 


mag 


1 


19.91 


19.24 


19.35 


5 


22.84 


21.40 


21.63 


2 


21.10 


20.23 


20.37 


6 


23.04 


22.46 


22.88 


3 


22.07 


20.47 


20.57 


7 


22.98 


20.43 


20.35 


4 


21.27 


20.71 


21.52 


8 


24.05 


23.32 


23.23 



errors are ±0^03 in the B and R c , and up to ±0^05 in the H a bands 



The B and R c magnitudes of the stars chosen as secondary 
standards were measured at the reference-night images cali- 
brated with the primary standards as described above. The stars 
are marked by numbers in Fig.Q] and their magnitudes with er- 
rors accounting for the zero-point errors shown in brackets of Eq 
(1) are listed in Table [7] 

Using the secondary standard magnitudes, the resulting zero- 
points obtained for the images summed over the whole 2003- 
2004 period are £=27.09+0.03 and fi f =27. 37+0.02. For the 
2001 period the resulting B band zero-point is 26.93+0.03. We 
also estimated 3<x detection limits of a point-like object in the 
co-added images. In the B band it is ~28™2 for the both 2001 
and 2003-2004 periods, and ~28™6 for the sum of the two peri- 
ods. It is ~27™9 in the R c band. There are no significant differ- 
ence between the B band formal detection limits for both peri- 
ods. However, an "eye visibility" of a 3<x object is better for the 
2001 period. This is due to better seeing conditions and a higher 
FORS1 resolution setup during this period (see Sect. [2J. 

7 http://www.eso.org 
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Stellar magnitudes Mj can be transformed into absolute 
fluxes Fj (in erg s crrT 2 Hz -1 ), whenever necessary, using 
standard equations 

log Fj = -0.4 (Mj + M°) , (5) 

with the zero-points provided by Fukugita (1995): M° B = 
48.490, M° R = 48.800. 

3.3. Flux calibration of the H image 

Flux calibration in the H a narrow band (/l=6563A, 
FWHM=60A) image was done using the spectrophoto- 
metric standard HZ44 observed the same nights as the pulsar. 
The flux Fff a and AB-magnitude M^ B of a source in this band 
are defined as 

/xl0 (Q.4^, eC z) erg 
b H a ~ L ; 5 TT~ (t, -> 

t exp cm 2 - s Hz 

Mj£= -2.5 log(F H J-4S.6, (7) 

where / is the measured flux of the source in counts and C the 
calibration constant. 

The calibration was performed using the night of 19 January 
2004 as a reference. From count-rate variations of stars in the 
pulsar field with the airmass, we derived a mean atmospheric 
extinction factor in the H a band kn n = 0.055 + 0.017. Then the 
calibration constant for this night was derived to be 7.66xl0~ 30 . 
The H„ magnitudes of the secondary standards measured using 
the reference night image are listed in Table [7] Based on these 
magnitudes the derived 3cr flux detection limit for a stellar object 
in the summed H^ image is l.lxl0~ 30 erg crrT 2 s _1 Hz _1 , and it 
is 5.5xl0~ 31 erg cirT 2 s _1 Hz -1 arcsec -2 for a surface brightness 
of an extended object. 

4. Results 

4. 1 . Searching for the pulsar counterpart in the B and R 
bands 

Zoomed fragments of the field containing the pulsar in the B 
and R bands for the 2001 and 2003-2004 periods are presented 
in Fig. |2] For comparison we also show the X-ray image of the 
same fragment obtained with the Chandra/ACIS-S in February 
20050 where a candidate X-ray pulsar counterpart has been 
found (Kargaltsev et al. I20061 >. The expected positions of the 
pulsar for each observing epoch, together with the 3cr position 
uncertainty ellipses from Tables [5] and [6] are marked with the 
pulsar proper motion direction. 

There are no significant objects within the pulsar position 
error ellipse in the 2001 B-band image. The nearest detected ob- 
ject, marked A, is outside the ellipse and lies 1'.' 1 away from the 
pulsar 2001 position in a direction not coinciding with the pulsar 
proper motion. Its B magnitude is 28.0 + 0.3 or only about 3.5<x 
of the formal detection limit. The outer contour of this object, 
where it merges with backgrounds, is overlaid on the 2003 B- 
band image. As seen, it only partially overlaps the pulsar error 
ellipse at the epoch of 2003, although its center is very near to 
the ellipse border. 

At the same time, within the pulsar position uncertainty el- 
lipse in the 2003 B-band image, we find another faint object (B) 

8 The X-ray data were retrieved from the Chandra archive (Obs. ID 
5519, 2005-02-23. 18ks exposure, PI G. Garmire. 



with a similar magnitude B = 28.1 ± 0.3. By a good positional 
coincidence with the expected pulsar coordinates, the object B 
can be considered as a candidate pulsar optical counterpart. The 
object A also can be resolved in this image within its contour 
of the 2001 epoch (see below). However, it sits in the east wing 
of the candidate and is not visible in the 2003 B-band image of 
Fig- El where the count range is chosen in a such way as to un- 
derline the presence of object B. The positions and magnitudes 
of objects A and B are listed in Table [8] 

A thorough inspection of the region containing objects A 
and B in the 2001 and 2003 B-band images was made with us- 
ing individual exposures. It confirmed that both objects are real 
but not artifacts caused by, e.g., a poor cosmic ray removing or 
flat-fielding. Owing to their closeness and faintness, these ob- 
jects may represent bright parts of the same unresolved extended 
background feature. To verify this possibility we considered the 
changes in spatial profiles of the region containing both objects 
from one observing epoch to another. Within positional uncer- 
tainties, the objects almost have the same declinations, while 
their right ascensions differ by about one arcsecond. This allows 
us to compare only ID E-W spatial profiles at a fixed declina- 
tion. 

In Fig. [3] we present the ID profiles extracted from the 2001 
and 2003 B-band images, shown in Fig. [2] along a horizontal slit 
with PA=90° and DEC=+ 15:51:1 1.49. The slit length and width 
are 4'.'6 and Of.' 6, respectively, and the centers of both objects 
A and B are within this slit. The coordinate origin of the hori- 
zontal axis in Fig. [3] corresponds to the eastern edge of the slit 
with RA=1 1:36:03.317. Only object A is seen in the 2001 profile 
(top panel). The 2003 profile is significantly wider, about twice, 
and has an asymmetric shape suggesting two sources of different 
intensities with overlapping profiles (bottom panel). Its peak is 
shifted towards west and corresponds to the position of object B, 
while object A is probably responsible for the "shoulder" in the 
east profile wing. 

The objects are too faint to be resolved reliably as two point- 
like sources with standard PSF modeling and subtraction tools. 
To analyze the situation, we applied a simplified approach based 
on the fact that only a single object A is resolved in the 2001 
B-band image and profile. We used its profile shape as a tem- 
plate of a single source profile to fit the suggested A + B blend at 
the 2003 period by a sum of the emissions from two sources, as- 
suming that B has the similar profile shape. Within uncertainties, 
this model is in a good agreement with the 2003 observations if 
the peak intensity of object A is by a factor of 1 .5 lower than 
that of B. This corresponds to about a 0.4 magnitude difference, 
which is within the uncertainties of the objects' magnitude mea- 
surements. The distance between A and B is about 078-l'.'0. This 
simplified approach supports the interpretation that A does not 
change its position and brightness significantly from the 2001 
epoch to 2003, while B appears in the 2003 image as an addi- 
tional source roughly of the same brightness, and by positional 
coincidence it can be associated with the pulsar. We discuss this 
and other possibilities of the A + B interpretation in Sect. [5] 

Neither of the objects A and B are resolved in the B-band 
image, which is obviously shallower due to a shorter integration 
time (cf. Sect 2). Nevertheless, the absence of a "red" source 
at the expected pulsar positions allows us to constrain the color 
index for the candidate detected in the B band as B - R < 0.5. 

To better compare the 2005 X-ray and the 2003 B-band im- 
ages we overlaid the contours of the optical image onto the X-ray 
image. The pixel scale of the ACIS-S image is about 0f.'5, which 
is comparable to the proper motion shift of the pulsar between 
the 2003 and 2005 epochs, but about three times smaller than 
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Fig. 2. Fragments of the PSR Bl 133 + 16 field images obtained in the optical B and R c bands with the VLT/FORS1 and in 0.2-10 
keV X-ray range with the Chandra/ ACIS-S at different epochs, as indicated in the plots. All images are smoothed with ~1" Gaussian 
kernel, and the fragment sizes are ~14"xl4". The + and x mark the expected radio positions of the pulsar at the observing epochs 
in the optical and X-ray images, respectively. The ellipses show 3cr position uncertainties, and the arrows show the direction of the 
pulsar proper motion and its path over about 20 yrs. The shift in the pulsar from 2001 to 2003 is indicated in the R c image. The 
contour of object A seen in the 2001 B-band image is overlaid on the 2003 B-band image. The contours on the X-ray image are from 
the 2003 B-band image. Note that object B in the 2003 fi-band image nicely fits the position of the candidate X-ray counterpart of 
the pulsar in the ACIS-S image. 



the FWHM of the ACIS PSF Within these uncertainties and a 
nominal ~1" Chandra/ ACIS pointing accuracy, the contour of 
the possible optical counterpart, object B, fits the position of the 
candidate X-ray counterpart nicely, which is the brightest object 
in the X-ray image. This implies, at least, that we probably see 
the same object in the optical and X-rays. The position of the 
X-ray candidate with 072 accuracy coincides with the expected 
radio position of the pulsar at the epoch of the Chandra observa- 
tions (Kargaltsev et al. 2006). We also note that counterparts of 
some other optical objects in the field may be found in the X-ray 
image and vice versa. 



4.2. Possible pulsar tail in H a and X-rays? 

The 2004 // ff -image of the same field and the difference in the 
H a and normalized /?-band images, where most of the contin- 
uum emission is subtracted, are shown in Fig. [4] The contours of 
the ACIS - S X-ray image, which is also shown here but with a 
slightly different gray scale than in Fig. [2] are overlaid on the 
(H a - /?)-image and vice versa. As in Fig. [2] all images are 
smoothed with the Gaussian kernel of about 1", which roughly 
corresponds to the seeing of the resulting optical images. 

In the H a image we find a faint point-like source (C), which 
lies within the pulsar position error ellipse. It is near the de- 
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Table 8. The parameters of the optical objects detected around the expected position of PSR Bl 133+16 



Object 


RA 


DEC 


epoch 


magnitude 


label 


hh mm ss 


dd mm ss 






A 


11 36 03.22+0.04 


+ 15 51 11.16+0.60 


2001 


B= 28.0 + 0.3 


B 


11 36 03.16+0.04 


+ 15 51 11.49+0.60 


2003 


B= 28.1 ±0.3, B-R < 0.5 


C 


11 36 03.15+0.04 


+ 15 51 10.86+0.60 


2004 


Hf = 26.3 ± 0.4 
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Fig. 3. The observed E-W spatial profiles of the regions contain- 
ing objects A and B extracted from the 2001 (top) and 2003 (bot- 
tom) B-band images of Fig. [2] along the slit with PA of 90° and 
Dec=+15:51:11.49. The slit length is 4'.'6 and its width is 076. 
The coordinate origin of the horizontal axis corresponds to the 
eastern edge of the slit with RA=1 1:36:03.317. The background 
was subtracted. Vertical dots mark the positions of A and B de- 
fined from the peaks of the profiles in the top and bottom pan- 
els, respectively. The short dashed line in the top panel shows 
a smoothed A profile that is used as a template to fit the A + B 
profile in the bottom panel as a sum of the two sources. The re- 
sulting fit is shown by the long dashed line, while short dashed 
lines show the contributions of A and B. 



tection limit and its star-like H^ B magnitude is 26.3+0.4 (2.7<x 
formal detection limit), which corresponds to the flux Fn a ~ 
1.05 x 10~ 30 ergs cirT 2 s -1 Hz -1 or 1.5 x 10" 6 photon crrT 2 s -1 . 
The object is not detected in the B and R c bands, while its possi- 
ble artifact origin in the narrow band was ruled out after a careful 
inspection of individual H a frames. 

We do not see any characteristic extended bow shock struc- 
ture around the pulsar, as is produced, for instance, around an- 
other high-velocity pulsar PSR B2224+65, known as the Guitar 
nebula (Cordes et al. I1993I I. Instead of that, we see a faint and 
rather clumpy H a emission most likely associated with the re- 
combination of a heated ambient matter around background ob- 
jects of the field. Some of them may have X-ray counterparts. 
However, there is one faint structure (about 4.5" length) south 
of object C. This apparent structure is roughly aligned with the 
pulsar proper motion direction and located behind its 2003 posi- 
tion. It is better seen in the (Hq.-R) image where it merges with 
object C. There are no signs of the structure in the broad band 



optical images. At the same time, in the X-ray image, we also 
see a marginal tail-like structure behind the pulsar counterpart 
candidate, though of smaller spatial extent. It was not mentioned 
by Kargaltsev et al. 2006 who only reported on the detection of 
the X-ray pulsar candidate. If both structures are not background 
fluctuations, they can be considered as cometary-like pulsar tail 
candidates. 



5. Discussion 

In Table [8] we collected the coordinates and magnitudes of the 
objects detected around the radio positions of the pulsar in each 
VLT observing epoch. Given that sources B and C are located 
within 3<t pulsar position error ellipces related to respective ob- 
serving epochs (cf., Tables|6]and[8]l, they can be associated with 
the pulsar. Object B, emitting in a broadband range, can be con- 
sidered as a candidate pulsar optical counterpart, while C can 
be a sign of the pulsar interaction with a clumpy ambient ISM 
seen via the H a recombination line. Object A is outside of the 
3cr-error ellipses for both observing epochs and we qualify it as 
a background source. If object B detected in the 2003 B-band 
image is the real counterpart of the pulsar, then why do we not 
resolve it in the expected place on the 2001 B image? 

One of the reasons is the faintness of the object, which is re- 
solved at a level very close to the estimated 3-cr detection limits 
(see Sect. 3.2). To estimate the probability of detecting an object 
with S/N ~ 3 in our data, we selected almost a homogeneous 
subset of 2001 B-band images, B12-B17, which provides a max- 
imal number of frames with the same seeing, background, and 
atmospheric extinction (see Table |3j. Then, we selected about 
100 faint objects chosen randomly in the sum of these 6 frames 
and checked that these objects are detected in the each of the in- 
dividual frames. That was done in two ways. In the first, we used 
the optimal aperture, with 3 pixel radius, and an object was ac- 
cepted as detected at an individual image if its magnitude error 
was less that 0.36, which corresponded to S /N ~ 3. In the sec- 
ond way we estimated the object by eye in the images smoothed 
with the Gaussian kernel. In the first case the probability of de- 
tecting an object with S/N ~ 3 was 60% - 70%. In the second, 
less formal case, it was about 80% - 90%. This means that, if 
we have two sets of VLT observations in the same band with ap- 
proximately equal conditions and exposure lengths and detect a 
faint object with S /N ~ 3 in one of them, the probability of los- 
ing this object in backgrounds at the second set of observations 
by pure count statistic is about 10% - 40%. Accounting for not 
quite homogeneous conditions for the 2001 and 2003 periods, 
a conservative estimate for losing the suggested pulsar optical 
counterpart at the expected position at the 2001 epoch is about 
50%. 

In addition, a slightly higher background, e.g., due to the 
closeness of a nearby bright object east of the pulsar path (cf. 
Fig. 2) or to a small increase in the interstellar extinction (e.g., 
by AE(B - V) ~ 0.04) for the line of sight towards the 2001 pul- 
sar position can drown the candidate in backgrounds in the 2001 
B image. 
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Fig. 4. The same as in Fig.[2]but for the H a and (H„-R) images, with the contours from the X-ray image overlaid. The X-ray image 
with the (H„-R) image contours overlaid is also shown for comparison. 



The pulsar optical brightness variation cannot be excluded 
also, since it is known to show a pulse nulling phenomenon in 
the radio range (e.g., Bhat et al. 2007 and references therein). It 
spends up to 15% of its time in the null state, and shows epochs 
when its radio flux changes by a factor of two. There are still 
no direct simultaneous radio and optical observations of pulsars. 
However, one cannot rule out a link between pulsar activities in 
these domains since both emissions are governed by complicated 
NS magnetosperic processes. The 2001 observing period of one 
week duration is much shorter than the 2003 one spanning 8 
months and a chance to meet a comparable flux depression phase 
that might completely hide the counterpart in the 2001 epoch is 
higher. 

We also cannot exclude the possibility that B is simply a time 
variable background object or a part of an extended unresolved 
variable feature that includes A and B as its relatively brighter re- 
gions. Another possibility is that A and B represent the same ob- 



ject displaced between 2001 and 2004 epochs by its own proper 
motion. However, our analysis has shown (Sect. 4.1) that A and 
B are likely to be two single objects and A does not show any sig- 
nificant variability or proper motion towards the position of B to 
explain the variability of the latter. We have also found no other 
variable objects in the pulsar neighborhood at the time base of 
three years. Possible variability of B can be checked only by fur- 
ther observations. If they will show that B survives and moves 
consistently with the pulsar this will be a strong proof that we 
see the real optical couterpart. 

A good positional coincidence of B with the candidate X-ray 
pulsar counterpart suggests that we see the same object in the 
optical and X-rays. Its color index, B - R < 0.5, is compatible 
with the indices of Vela pulsar, PSR B0656+14, Geminga, and 
PSR B0950+08, which all have B - in a narr ow ran ge of 0.46 
- 0.6 (Mignani & Caraveo l200~fl Shibanov et al. 2006; Zharikov 
et al. 11004'). 
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Fig. 5. The evolution of the pulsar radio, optical and X-ray lu- 
minosities, and respective efficiencies with the characteristic age 
t. The top panel shows corresponding spin-down luminosities. 



Let us assume that the optical object B and the respective X- 
ray source are indeed the pulsar PSR Bl 133+16. The spectral fit 
of the Chandra/ ACIS-S data of the X-ray counterpart candidate 
by an absorbed power law (Kargaltsev et al. 2006) yields the 
non-thermal X-ray luminosity of the pulsar Lx = 5. 01 (+4. 30 / - 
2.41) x 10 28 in the 2-10 keV range at the distance of 357 pc, 
or logLx = 28.70(28). The corresponding X-ray efficiency is 
log77x = -3.24(28), which is not exceptional and lies near the 
efficiency range of the two other old pulsars, PSR B1929+10 and 
PSR B0 950+0 8, detected in the optical and X-rays (Zharikov et 
al. 120041 [2006b . 

Using i]x and the empirical relation between the non-thermal 
optical and X-ray efficiencies of ordinary pulsars reported by 
Zharikov et al. (2004, 2006), we would expect to find the 
PSR Bl 133+16 optical counterpart with a magnitude in a range 
of about 29-30 and with the optical efficiency log r\opt ~ -6.2. 
The proposed candidate is a magnitude brighter and the respec- 
tive optical luminosity and efficiency, logL^ = 26.76(17) and 
log rjopt = -5.18(17), are also about one order of magnitude 
higher than is expected from the reported relation. However, this 
may simply reflect inaccuracies in the empirical relation, which 
is based on a very limited sample of the optically detected pul- 
sars. 

To study how this discrepancy can modify the relation and 
the evolution tendencies considered by Zharikov et al. (2004, 
2006), we included the proposed candidate in the sample. The 
results are shown in Figs. [5] and [6] As seen (Fig|5]l, our candi- 
date does not invalidate the previous conclusion regarding the 
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Fig. 6. Relations between the optical j]o pt and X-ray r\x efficien- 
cies in the B-band and 2-10 keV ranges. Solid line shows the 
best linear regression fit for a sample of pulsars from Zharikov et 
al. (|2006). The dot-dashed line is the similar fit only for the set of 
the pulsars that includes PSR Bl 133+16 counterpart candidate 
and excludes the youngest Crab and PSR J0540-69. Numbers 
are the line slopes with their errors in brackets. 



high optical and X-ray efficiencies of old pulsars, which be- 
come comparable to those of the young Crab-like pulsars. At 
the same time, its position on the rjopi-rjx plane significantly 
deviates from the previous relation (Fig. [6j. Two of its neigh- 
bors, PSR B1929+10 and PSR B0950+08, also show apparent 
deviations. However, if we exclude from the sample the two 
youngest pulsars, the Crab pulsar and PSR B0540-69, we ob- 
tain a new relation that fits the rest pulsars nicely, including the 
PSR Bl 133+16 candidate. This is not a surprise, since noticed 
Zharikov et al. (2004) that the young pulsars are distinct from 
the older ones by their X-ray and optical efficiencies and dis- 
placement on the T]opt-E and rjx-E planes. Therefore, in both 
respects the proposed candidate is not outstanding. That can be 
considered as additional evidence that it is the real counterpart 
ofPSRB1133+16. 

Despite the high transverse velocity of PSR Bl 133 + 16, we 
did not detect any developed Balmer bow shock nebula of the 
pulsar. This suggests low ISM density in the pulsar neighbor- 
hood that is compatible with its high galactic latitude «69", low 
E(B - V)=0.04 towards the pulsar, and a small hydrogen col- 
umn density A^=(l-2)xl0 20 cnT 2 (Kargaltsev et al. 2006]! es- 
timated from the measured pulsar dispersion measure of 4.86 
cirT 3 pc and the spectral fit of the candidate X-ray counterpart. 
Nevertheless, the emission detected at about the 3<r detection 
limit in the pulsar position error ellipse in the Hq, image (object 
C) may come out from the brightest part of the bow shock head 
and/or from a clump of the ISM with a higher density. The low 
S/N ratio does not allow us to conclude this confidently. The ab- 
sence of object C in our B and R broad band images prabably 
excludes its background origin. Accounting for deep B and R 
magnitude upper limits and a high galactic latitude with low 
E(B - V), object C certainly cannot be a star from our Galaxy. 
Any background extragalactic object is also very unlikely since 
it would be detected first of all in our deep broad band images. 

Instead of the bow shock, we note two marginal but likely 
spatially correlated extended features in the Hq, emission and 
X-rays, which are reminiscent of a "tail" behind the pulsar. If 
the X-ray and the H„ tails are indeed associated with the pul- 
sar, this would present the first example of when the pulsar tail 
is detected simultaneously in the H a and X-rays. The combina- 
tions of the X-ray tails and H a bow shocks for a few pulsars 
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have been reported in the literature (see, e.g., Zavlin 120061 for 
a short review). Most recently such a tail has been reliably de- 
tected linX-rays with the XMM behind PSR B 1929+ 10 (Becker 
et al. |20061 l, which is also a member of the sample considered 
above. However, despite these efforts, no counterparts of the 
tails have been detected so far in the optical range to our knowl- 
edge. The formation of the H a tails does not look unreasonable 
when a compressed pulsar wind behind a supersonically moving 
NS cools via X-ray emission that ionizes the ambient matter on 
larger spatial scales, which then recombines producing H a pho- 
tons. An apparent offset of the X-ray tail candidate from the Hq, 
elongated structure in our case (cf. the images with contours in 
the H ff -R and X-ray panels in Figj4j is consistent with such a 
scenario. 

The suggested candidate optical and X-ray counterparts of 
PSR Bl 133+16 and its tail can be easily verified by a followup 
imaging of the field in the optical and X-rays on the basis of 
a few years. Given the pulsar proper motion, the candidates, if 
they are the real counterparts, will be shifted from the detected 
positions roughly to the north by about 1'.'4-1'.'8 by the end of 
2008. Such a shift can be reliably measured using the subarc- 
second spatial resolutions of 8m-class optical telescopes and the 
Chandra X-ray observatory. 
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